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ABSTRACT 

We observed the energetic binary Cygnus X-3 in both quiescent and flaring states 
between 4 and 16 /j,m using the ISO satellite. We find that the quiescent source shows 
the thermal free-free spectrum typical of a hot, fast stellar wind, such as from a 
massive helium star. The quiescent mass-loss rate due to a spherically symmetric, 
non-accelerating wind is found to be in the range (0.4-2.9) x 10~ 4 Mq yr -1 , consistent 
with other infrared and radio observations, but considerably larger than the 10~ 5 
Mq yr _1 deduced from both the orbital change and the X-ray column density. There 
is rapid, large amplitude flaring at 4.5 and 11.5 /zm at the same time as enhanced 
radio and X-ray activity, with the infrared spectrum apparently becoming flatter in 
the flaring state. We believe non-thermal processes are operating, perhaps along with 
enhanced thermal emission. 

Key words: stars: individual: CygX-3 - stars: Wolf-Rayet - binaries: close - stars: 
mass-loss - infrared: stars - X-rays: stars 



1 INTRODUCTION 

The Infrared Space Observatory ISO (Kessler et al. 1996) 
with its on-board complement of cameras and spectrome- 
ters, provided a rare opportunity to study sources at wave- 
lengths not accessible from the ground. We present photome- 
try and images of the Cygnus X-3 field made using ISOCAM 
(Cesarsky et al. 1996), the camera on ISO, in the wavelength 
range 4 to 16 /im. We report the first spectrum of this source 
in this wavelength range. 

Cygnus X-3 is an enigmatic and unusual X-ray binary 
which was discovered >30 years ago (Giacconi et al. 1967). 
Observations of it have been made from radio to gamma rays 
and in every wavelength bands the system exhibits unusual 
behaviour. In X-rays and infrared, the emission is partially 
modulated by a 4.8 h period (Mason, Cordova & White 
1986), and this is believed to be the orbital period of the 
system. 

In 1991, van Kerkwijk et al. observed Cyg X-3 in the 
infrared with the CGS4 spectrometer on UKIRT and discov- 
ered He I and He n together with N emission lines which were 
modulated at the 4.8 h period (van Kerkwijk et al. 1992; van 
Kerkwijk et al. 1996). They concluded that the lines origi- 
nated from the companion star, either a Wolf-Rayet of the 
WN7 class, or an earlier WN4-6 in an outburst phase. This 
was investigated in more detail by Fender, Hanson & Pooley 
(1999, hereafter FHP99), who took spectra during a variety 



of activity states and confirmed that the underlying quies- 
cent spectral type is of an earlier WN4-5 Wolf-Rayet star. 

Galactic Wolf-Rayet stars have high velocity winds, up 
to 2500 km s" 1 (Crowther, Hillier & Smith 1995; Bohannan 
& Crowther 1999). The Wolf-Rayet wind from Cyg X-3 is 
believed to be the main contributor in providing mass and 
angular momentum loss which increase the orbital period 
(P = -6.6 x 10~ 10 s s"\ Kitamoto et al. 1995). This pe- 
riod increase is consistent with a mass- loss rate of 10 _6 Mt 
Mq yr _1 , where Mt is the total mass in the system (Kita- 
moto et al. 1995; Matz 1997). Kitamoto et al. (1994) mod- 
elled changes in the radio activity, infrared and X-ray states 
caused by changes in the density of the WR wind (Table Q). 
The multi- wavelength predictions of Kitamoto et al. (1994) 
were tested by FHP99, who find that during an outburst 
the infrared K band spectrum becomes dominated by strong 
twin-peaked He I emission lines. These lines show strong 
V/R variability correlated with the orbital phase. Due to 
the small size of the binary (~ 5 Rq separation, 4.8 hour 
period) the observed line flux could not arise from a region 
of this size as the environment would be too hot. FHP99 
conclude that the He I emission lines are emitted from a re- 
gion which is significantly larger than the binary separation, 
and that the variability of these lines indicate an asymmetric 
emission region, possibly from a flattened disc-like wind. 

The compact object in Cyg X-3 is less well defined. 
Cherepashchuk & Moffat (1994) suggest that the high X-ray 
luminosity (greater than the Eddington luminosity for a 1 
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Table 1. Radio, infrared and X-ray variability based on a chang- 
ing density wind, from Kitamoto et al. (1994). The term 'mod- 
ulation' refers to the continuum. The states marked with an as- 
terisk (*) have been expanded on and confirmed by Fender et al. 
(FHP99). 
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Mq object) implies that the compact object is a black hole. 
Schmutz, Geballe & Schild (1996) related time-variations 
of infrared spectral lines to an orbital velocity and derived 
a mass for the compact object in the range of 7-40 Mq, 
assuming realistic values for the mass of the Wolf-Rayet 
star and the inclination of the system. However, their model 
does not account for all the relative phasing of the X-rays 
and infrared lines. Ergma & Yungelson (1998) point out that 
if the progenitor to the binary is two massive helium stars 
then the mass-loss and common envelope phase of these two 
stars prevents the formation of a black hole & Wolf-Rayet 
binary with an orbital period less than several days. This 
ambiguity does not resolve itself if the compact object is 
a neutron star; Ergma & Yungelson (1998) also comment 
that the propeller motion of a spinning neutron star would 
prevent wind accretion in the system, thus removing the 
jet-creation mechanism. 

With a high mass, high mass-loss (wind) object in close 
proximity to a compact object, a high degree of interaction 
takes place. Radio observations at wavelengths ranging from 
1-20 GHz have shown Cyg X-3 to undergo many states of 
emission and are used to define the state of Cyg X-3. For 
most of its time, the system occupies a quiescent radio state, 
with 8.3 GHz flux around 100 mjy (Waltman et al. 1994). 
Major flares (8.3 GHz flux > 1 Jy) occur at random in- 
tervals every 18 months or so and relaxation after a flare 
can last up to several weeks (Waltman et al. 1995). More 
common are minor flares (8.3 GHz flux 200-1000 mjy) and 
these usually occur during a flaring period typically several 
weeks duration, with several small flares occurring in close 
proximity to each other (Waltman et al. 1996). Minor flares 
can reach peak emission in a few hours and decay over sev- 
eral hours (Ogley 1998). The minor and major flares have 
been shown to be intrinsically different in their production, 
rise and relaxation times (Newell 1996). Preceding a ma- 
jor flare, the system enters a quenching state with fluxes at 
8.3 GHz dropping to around 10 mjy. The system returns to 
quiescence then flares (Waltman et al. 1996). 

The main mechanism for radio flares is thought to be 
associated with the acceleration of relativistic electrons in 
the base of jets. Synchrotron emitting plasmons have been 
detected, travelling at ~ 5 mas d^ 1 on both arcsec and mas 
scales (Spencer et al. 1986; Mioduszewski et al. 1998). There 
is no confident determination of the velocity of ejecta: while 
a velocity of 0.3c is generally taken as the speed on the 
larger arcsec scale (e.g. Spencer et al. 1986), other highly- 



Table 2. ISOCAM filters and integration times used in the qui- 
escent observations. 



Filter 


Wavelengths 


Integration 


name 


range (fJ.va) 


time (s) 


LW 1 


4.10-4.90 


161 


LW 4 


5.50-6.50 


161 


LW 5 


6.50-7.00 


119 


LW 6 


7.00-8.50 


144 


LW 10 


8.60-14.4 


321 


LW 3 


12.1-16.9 


311 



relativistic interpretations are possible on the mas scale 
(e.g. Ogley, Bell Burnell & Newell 1997; Newell, Garrett & 
Spencer 1998). 

The Wolf-Rayet contributes not only to the enigmatic 
qualities of the system, but also a significant and varying 
opacity. The true nature of the secondary is a top priority 
for the understanding of the system. Observations at the 
ISO wavelengths provide a unique opportunity to study the 
Wolf-Rayet wind, and spectra obtained can track its true 
nature. 



2 ISOCAM OBSERVATIONS 
2.1 In quiescence 

Quiescent observations were taken on 1996 April 07 (MJD 
50180) with 6 long-wavelength (LW) filters on the CAM in- 
strument. Cyg X-3 is in the galactic plane at a distance of 
around 10-kpc and so is heavily obscured by dust along our 
line of sight. Around the object, emission from cirrus is par- 
ticularly strong at longer wavelengths so our observations 
were focused in the 1-16 um range. The wavelengths chosen 
are presented in Table |2j We observed for a total of 2895 
s in 6 filters in order to obtain a rough spectrum of the 
source. A 6 arcsec pixel size was used for the observations 
and single frames, stare-mode, were taken with the 4 shorter 
wavelength filters, whereas the 12 and 15 /im images were 
taken as 3 x 3 mosaics. The integration time for individual 
frames was 2.1 s and the total integration times for all filters 
are given in Table ^ 

Analysis of the ISO data was performed at the ISO 
data centre in Vilspa, Spain. Each image had a number of 
corrections applied to it and these included: (a) dark cor- 
rection, (b) de-glitching the data to remove transients, (c) 
stabilisation of pixels due to memory effects, (d) a flat-field 
correction, (e) jitter correction, (f) combination of images 
into a stare image or raster map, (g) field of view distortion 
and finally (h) conversion from engineering units of adu g~ 
s _1 to Jy. 

The corrections that were applied to the data either 
used standard library frames, or software routines developed 
by the ISOCAM Consortium. The processes which used li- 
brary fields were the dark correction (van Buren & Kong 
1996) and the flat-field correction of stare images. The other 
corrections used standard software routines and included a 
multi-resolution median transform for the de-glitching pro- 
cess; a 90 per cent stabilisation method to correct for mem- 
ory effects; standard field of view distortion; and standard 
conversion of adu g _1 s _1 to Jy with code formulated on 
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1998 April 05. No jitter correction as the average jitter re- 
ported of 0.5 arcsec was much less than our pixel field of 
view of 6 arcsec. For a detailed discussion see Ogley (1998). 

The reduction of mosaic images was essentially identical 
to stare images, but with the following exceptions. The flat- 
field correction did not use standard frames, but built up 
a flat-field using the background pixels in the mosaic. The 
flat-fielding using standard frames is the least robust of the 
data analysis stages, and can lead to significant increases or 
decreases in flux at the edges of an image. Memory effects 
due to pixel stabilisation can create a cross-shaped effect on 
mosaic images. This was reduced by sequentially masking 
the residual pixels as the mosaic is built up. 

2.2 In outburst 

Three epochs of ISO data were acquired to observe Cyg X-3 
during a major flare on 1997 Jun 15, 18, 21; MJD 50613.7, 
50615.4, 50619.9. Fig. [j] shows photometry at 8.3 GHz from 
the Green Bank Interferometer (GBI) with the times of our 
ISO observations in both quiescent and flaring states shown 
by bullet points, one epoch during quiescence, and three 
epochs during outburst. Each flaring epoch consisted of a 
stare observation with CAM using the LW1 filter (4.5 /im), 
a mosaic observation with CAM using the LW10 filter (11.5 
/im) and a PHOT observation. Due to instrumental difficul- 
ties only the CAM data is presented. The procedure used in 
reducing the data was identical to the quiescent epoch. 



3 RESULTS 
3.1 The Images 

3.1.1 Digitized sky survey 

There is a striking difference between the ISO images and 
an optical digitized sky survey (DSS) image of the same 
region. The DSS image, Fig. ^, shows many more stars than 
are visible on the ISO images (Fig. |^). Cyg X-3 is undetected 
on the DSS plate as it has Ay ~ 30 mag of extinction, and 
an estimated V band magnitude of V — 29 ± 1 (Wagner 
et al. 1990). The radio position of the source is represented 
by a cross, at the centre of the image. Stars marked KMJ 
1315, 1327, 1328 (Kobulnicky, Molnar & Jones 1994) and 
star Westphal 'A' (Westphal et al. 1972) are all identified. 
KMJ 1328 is Westphal's star 'C with V = 15.03. One can 
see that, as is expected along the galactic plane, there is an 
abundance of sources at various magnitudes. 

3.1.2 ISOCAM-quiescence 

In contrast to the DSS image, Fig. ^ shows contour images 
at the wavelengths 4.5-14.5 pm. The images have been flat- 
fielded and dark-subtracted using the techniques described 
above. The data were exported and imaged using the NRAO 
Astronomical Image Processing System (AIPS), and a rota- 
tion of +99? 17 was applied to each image. 

The absolute position of the ISOCAM images is un- 
known due to mis-alignment after a filter change. The co- 
ordinates of a source can be offset from the true position 
by as much as 2 pixels, or for our observations, 12 arcsec. 



GREY: DSS image 
CONT: ISO at 4.5 urn 
2 4 6 8 10 12 




20 32 35 30 25 20 15 

RIGHT ASCENSION (J2000) 
Grey scale flux ranges 1.7 13.3 Kilo 
Conl peak flux = 2.0591 E+01 mjy 
Levs= 1.0000E+00 ' ( 10.00,12.00,14.00, 
16.00, 18.00) 

Figure 2. Digitized sky survey image of the Cyg X-3 area. The 
box shows the orientation and position of the stare-mode ISO- 
CAM images. For further details see the text. 

Positions can be calculated from phase-referenced radio im- 
ages, optical sources or catalogued sources. The positions 
were calculated from a radio observation of Cyg X-3 using 
MERLIN at 22 GHz with a 30 mas beam size (Ogley et al. 
2000). The radio position has a Gaussian-fitted position of 
RA(2000) = 20 h 32 m 25*712, Dec. (2000) = +40° 57' 28'.'21, 
and the ISO positions of all sources are correlated to this. 

In some images (6.0, 6.75 and 7.75 fim in Fig. ^) there is 
a significant misalignment between the flat-field frame and 
the array which produces a gradient at the edge of the im- 
age. Other techniques to flat field the data produced larger 
errors. As none of the sources of interest are near to the af- 
fected region, the photometry presented below is still valid. 

3.2 Photometry 

3.2.1 Quiescence 

Photometry on all four sources was performed using private 
code by Dr Siebenmorgen at Vilspa. An aperture of 12 arcsec 
was used on all sources and the flux densities are given in 
Table ^[ Spectra of the four sources are shown in Fig. ^. No 
data have been de-reddened. 

3.2.2 ISOCAM-outburst 

Photometry of the sources in the three flaring epochs were 
obtained in the same way, as in the quiescent observations, 
with a 12 arcsec aperture. The results for the photometry is 
presented in Table ^, for all sources. Cygnus X-3 is clearly 
variable compared to the field stars. 
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Figure 1. Radio photometry at 8.3 GHz from the Green Bank Interferometer. Displayed by the bullet points are the times of our four 
ISO observations, one epoch during quiescence, and three epochs during outburst. No data was available between MJD 50200 and 50400 
as the GBI was off-line for financial reasons. 



Table 3. Photometric fluxes of all four sources, during a Cyg X-3 quiescent state. A 12-arcsec aperture was used, corrected for the 
different point-spread functions of the filters. No data have been de-reddened. 



Star 4.5 fim flux 6.0 fim flux 6.75 /im flux 7.75 fim flux 11.5 fim flux 14.5 fim flux 

(mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 



KMJ 1327 


35.4 


± 


2.2 


50.2 


± 


3.6 


30.9 


± 


5.1 


30.6 


± 


4.9 


8.9 ± 


3.0 


6.5 ± 


2.7 


KMJ 1328 


26.0 


± 


1.8 


28.0 


± 


3.6 


22.0 


± 


5.7 


17.6 


± 


6.8 


8.0 ± 


2.6 


3.9 ± 


2.0 


Westphal et al. 'A' 


21.0 


± 


2.2 


22.2 


± 


3.2 


14.2 


± 


4.9 


22.4 


± 


4.9 


7.0 ± 


1.6 


4.3 ± 


1.5 


Cyg X-3 


18.2 


± 


2.5 


31.1 


± 


4.1 


20.6 


± 


4.9 


24.3 


± 


3.9 


15. 2± 


1.6 


18. 7± 


1.8 



Table 4. Infrared fluxes for the four sources against epoch during 
a Cyg X-3 flaring state. No source has been de-reddened. The top 
panel shows data at 4.5 (im and the bottom panel shows data 
at 11.5 (im. The orbital phase of Cyg X-3 is shown using the 
ephemeris by Matz (1997). 



4.5 fim flux (mjy) 
MJD 50613 50616 50619 

Date 1997 Jun 15 18 21 



Cyg X-3 phase 


0.85 


0.64 


0.54 


Cyg X-3 
KMJ 1327 
KMJ 1328 
Westphal et al. 'A' 


28.1 ± 2.3 
27.7 ± 1.9 
23.9 ± 2.2 
20.5 ± 2.4 


53.8 ± 3.1 

32.9 ± 2.4 
24.3 ± 2.0 
19.5 ± 2.3 


68.5 ± 3.2 

32.6 ± 2.5 
23.1 ± 1.9 
19.5 ± 2.6 


Cyg X-3 


11 

10.6 ± 0.6 


5 fim flux (m 
78.2 ± 0.8 


Jy) 

40.4 ± 0.8 



3.3 Spectra 

The Cyg X-3 data have been de-reddened using an extinc- 
tion coefficient in the J band (1.25-^m) of Aj = 6 (Fender 
et al. 1996), and an extinction law of A" 1 ' 7 (Matthis 1990). 



3.3.1 Quiescent spectrum 

The de-reddened data for Cyg X-3 in quiescence are plotted 
in Fig. | along with UKIRT data from 1994 July 17 taken at 
orbital minimum (Fender et al. 1996; Fender, Bell Burnell & 
Pooley 2000) . During the 1994 epoch observations a number 
of short-timescale flares occurred, but we have been careful 
only to use data at a time when the source was in quiescence. 
The spectrum shows a flattening at the longer wavelengths, 
with a spectral index of 3.0 between 1.25-2.2 fim and 0.4 
between 4.5-14.5 fim. 



3.3.2 Flaring spectra 

Table |B| shows the de-reddened photometry and spectral in- 
dices for the three observations. These spectra are also plot- 
ted in Fig. |. The UKIRT data from 1984 is from Fender 
et al. (1996) and is shown for two nights on 1984 Aug 07, 
09. Figure shows the flare spectrum with the quiescent 
spectrum subtracted. 



4 DISCUSSION 
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Figure 3. ISOCAM images of the Cyg X-3 field in quiescence. Images from top-left to bottom-right are at 4.5, 6.0, 6.75, 7.75, If. 5 and 
f4.5 /im. Photometry results are presented in Table | 
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Figure 4. Spectrum of the four sources, KMJ 1328, 1327, Westphal et al. 'A' and Cyg X-3. No data have been de-reddened. 



Table 5. Dc-rcddcncd fluxes for Cyg X-3 over the three flar- 
ing epochs. Data have been de-reddened using an extinction of 
A j = 6 (Fender et al. 1996), an extinction law which follows 
X~ l r (Mathis 1990) and the spectral index uses the definition 
S<xu a . 



Epoch 


4.5 fim flux (mJy) 


11.5 /im flux (mJy) 


a 


50613 


52.6 ± 4.3 


12.0 ± 0.7 


1.6 


50616 


101 ± 6 


87.8 ± 0.9 


0.1 


50619 


128 ± 6 


45.9 ± 0.9 


1.1 



4.1 The mass- loss rate 

In the ISO regime the spectrum becomes like that pre- 
dicted by Wright & Barlow (1975; hereafter WB75) for 
a non-accelerated, spherical wind. At long infrared wave- 
lengths, WB75 showed that free-free emission in a wind 
would display an S oc u 2 ^ z spectrum, while at longer radio 
wavelengths this is modified to produce a S oc i/ ' 6 spec- 
trum. Our de-reddened quiescent ISO fluxes can be fitted 
with a function S(v) = So (i^ /10 9 ) with the parameters 
So = 0.04 ± 0.06 and a = 0.61 ± 0.13 to give a reduced x 2 
value of 1.74. The So parameter is the flux at 1 GHz and its 
value around shows there is negligible expected free-free 
emission at radio frequencies. 

There is no evidence for infrared emission from dust, 
consistent with the absence of dust in WNs demonstrated 
by Cohen (1995). This means that the dust responsible for 
the soft X-ray halo emission reported by Predehl & Schmitt 



(1995), at an angular separation of at least 100 arcsec, can- 
not lie close to Cyg X-3. 

Using the WB75 model to obtain first-order wind pa- 
rameters, the free- free emission flux is given by 

2.3 XlO 10 ( M V /3 ( 2 s2/3 

S = — (vGHz-ycgZ ) mJy, (1) 

where Dk pc is the distance to the source in kpc, M is the 
mass lost in the wind in M© yr _1 , fi is the mean atomic 
weight per nucleon, Vao is the velocity of the wind at infinity 
measured in km s _1 , ^ghz is the observing frequency in GHz, 
7o is the number of free electrons per nucleon, g is the Gaunt 
factor and Z is the mean ionic charge. From our ISO data, 
the quiescent flux at a wavelength of 6.75pm (44,400 GHz) 
is 30 ± 10 mJy. Substituting these values into equation (1) 
together with a distance of 10 kpc, and assuming a value 
for the Gaunt factor of unity, gives an expression for the 
mass-loss rate of 

M = 3.2 x 10~ 8 ,Woc (y c Z 2 )~ 1/2 M yr" 1 . (2) 

The velocity of the wind in Cyg X-3 has been investigated 
by a number of authors, van Kerkwijk et al. (1992; 1996) 
observed Cyg X-3 using a high resolution spectrograph on 
UKIRT. Their observations were scheduled over a range of 
orbital phases and from this they were able to obtain a wind 
velocity of Voo — 1600 ± 300 km s _1 . More recently, the 
flattened wind model of FHP99 gives an outflow velocity of 
Woo ~ 1500 km s _1 . It is this latter value that we will adopt 
in this paper. 
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Figure 5. Quiescent and flaring spectra for Cyg X-3, observed by UKIRT at 1-4 (im and ISO at 4-15 fim. Included is UKIRT photometry 
from 1984 (Fender et al. 1996), and from 1994 at orbital minimum, from Fender et al. (2000). All data have been de-reddened using an 
absorption co-efficient of Aj = 6 (Fender et al. 1996) and a A~ 1,7 dependence (Matthis 1990). 



The other three variables which are unknown - fj,, 7 C 
and Z - are all related to the wind's composition. 

(a) A WN-type wind. Using a wind composition typical 
of late WR stars, /j, — 1.5, "/ c — 1, Z — 1 (Waltman et al. 
1996) and the flux at 6.75/im (44,400 GHz) the mass- loss 
rate is M a = 7.8 x 10" 5 M Q yr" 1 . 

(b) A partially ionised wind. The composition of this 
wind was used by Waltman et al. (1996) to explain the time 
difference between radio flares. This wind has the parame- 
ters of jj, = 4, 7 C = 0.5, Z = 0.5. The mass-loss rate using 
our infrared fluxes and these values is Mb = 2.9 x 10 -4 Mq 

yr" 1 . 

(c) A fully ionised wind, van Kerkwijk et al. (1996) used 
this type of wind when interpreting UKIRT data for the ob- 
servations of the Wolf-Rayet-like lines. This wind has the 
parameters of p = 4, -y c = 2, Z — 2 and using our infrared 
fluxes a mass-loss rate of M c = 3.6 x 10~ 5 Mq yr -1 is cal- 
culated. 

We note that Crowther et al. (1995) have found averages 
of fj, = 3.23, 7o = 1.010, Z = 1.015 for Galactic WN(L) 
stars (WN7-8) which would infer a mass-loss rate of M = 
8.1 x 10" 5 M Q yr -1 . For a WN(E) star (WN4-6) one might 
expect ^ ~ 4, Z = 1, 7 = 1 (e.g. Willis 1991) giving M ~ 
2.1 x 10" 4 M Q yr" 1 . 

Thus we find that for a single temperature, non- 
accelerating, spherical wind and for a variety of ionisation 
states, the free-free emission causing the infrared fluxes in 



the ISO wavelengths is due to a wind with a mass loss rate 
of (0.4-2.9) xlO" 4 M yr -1 . 

This value for the mass loss rate is towards the up- 
per end of the range observed for (single) WR stars, and is 
slightly higher than that predicted by Waltman et al. (1996) 
to be responsible for the delay between radio emission at a 
range of frequencies. However, it is slightly lower than is pre- 
dicted by van Kerkwijk et al. (1996) for their interpretation 
of the short wavelength infrared fluxes. Results from both 
papers together with our fluxes are given in Table ^. 

More significantly, there is a discrepancy in the value of 
M derived from the orbital period change and the value from 
our ISO results. Using the relationship P/2P = M/Mr and 
a value for P of -6.6 x 10" 10 s s" 1 , Kitamoto et al. (1995) 
calculated a value of M = 0.6 x 10" 6 Mr M yr -1 where 
Mt is the total system mass assumed to be of the order 
10 Mq. This is an order of magnitude lower than that from 
our ISO results (10~ 4 Mq yr -1 ). This difference between 
M derived from the orbital period change, and M derived 
from thermal emission observations (with the latter being 
the larger) is not unique to Cyg X-3. The mass-loss rate for 
the binary V444 Cyg (WR139) has been calculated to be 
6 x 10" 6 and 2.4 x 10~ 5 M yr" 1 by the P orb method and 
thermal emission method respectively (Underhill et al. 1990; 
Prinja et al. 1990). 
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Figure 6. The flare spectrum with the quiescent spectrum (as shown in Fig. a) subtracted. 



Table 6. The mass loss rate for the Wolf-Rayet star given by the 
various models discussed in the text. 



Method 


Ma 


Mb 






xlO 


" 5 M yr 


-l 


Radio photospheres 1 


< 0.25 


< 1.65 


< 2.7 


This paper 


7.8 ± 2.1 


29 ±8 


3.6 ± 1.0 


Near-IR fluxes 2 






> 12 


Orbital period change 3 


0.06 







1 Waltman et al. (1996) 

2 van Kerkwijk et al. (1996) 

3 Kitamoto et al. (1995) 



4-1.1 M from S - an examination of assumptions 

In obtaining M from the ISO data we have assumed that 
the wind is isothermal, at a constant velocity, shows spher- 
ical symmetry, stationarity and homogeneity. We have fur- 
ther assumed that the free-free flux is not seriously contam- 
inated by stars Z and D identified on the Cyg X-3 finder 
chart and within the ISO pixel size (Fender & Bell Burnell 
1996) nor by synchrotron emission from the tail of the radio- 
mm spectrum. The agreement with the Wright and Barlow 
model suggests that these assumptions (with the possible 
exception of spherical symmetry) are reasonable. We have 
assumed a distance of 10 kpc to the source; since M scales as 
Z} 3 / 2 , to effect a 10- fold reduction in the mass- loss rate, the 



distance would have to reduced to ~ 2 kpc. The extinction 
to Cyg X-3 is uncertain with 4.5 < Aj < 7.5. This corre- 
sponds to range of extinction at 6.75 yum of 0.26-0.43, and 
an uncertainty of only 10 per cent in the value of the flux. 

We have also assumed that the Gaunt factor, g, is unity. 
In fact it depends on Z, T and the observing frequency, and 
occurs in the expression for M as <?~ • The equation for the 
Gaunt factor given by Spitzer (1962), and subsequently used 
by Leitherer & Robert (1991) is valid for v p « v « kT/h, 
where v v is the plasma frequency. This condition is satisfied 
over the range of parameters we use here, and taking ex- 
treme values for the temperature of 10,000 to 350,000 K. 
Spitzer's equation gives values of g ranging from 3.4 to ef- 
fectively zero (at the lowest temperature and highest Z). 
Alternatively Sutherland (1998) determines the Gaunt fac- 
tor for the spectral range submillimetre to hard X-ray and 
for temperatures from 10 to 10 9 K. We believe that in his 
table 2 the labelling of the rows and columns has been in- 
terchanged, and that figures 2b and 2c have been switched. 
Sutherland's work then gives Gaunt factors that range from 
1.5 to 3.2. Sutherland considers photoionisation (which may 
be more appropriate to Cyg X-3) as well as collisional ion- 
isation and shows that the photoionisation figures for the 
Gaunt factor are typically a factor of 2 lower. Our assump- 
tion that the Gaunt factor is unity well represents the middle 
of the range obtained, but it is uncertain by a factor of ~ 3. 

Nugis, Crowther & Willis (1998) model the mass- loss 
rate in a wind that is clumped, and in which both the ion- 
isation structure and the clumping changes with radial dis- 
tance. Selecting values for their free parameters that best fit 
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for each star the flattening spectral slope from infrared to 
radio, they obtain mass-loss rates for clumped winds for 15 
WN stars. There is considerable scatter in their results, but 
they obtain an average of log M(clumpy) — log M(smooth) = 
— 0.2±0.3 (If). In other words, their clumpy wind model on 
average reduces the value of M to 0.6ig'j} of the M(smooth) 
value. Assuming Cygnus X-3 is close to average, then again 
this modification alone is insufficient to explain the discrep- 
ancy. 

FHP99 show that the most likely interpretation of their 
infrared spectra and their observation of an orbital modula- 
tion of V/R involves a non-spherical wind, with a flattened, 
disc-like enhanced wind, probably in the plane of the binary. 
They note that this also reconciles the high infrared emissiv- 
ity with an optical depth to the X-ray source that remains 
small (as long as the system is not viewed edge-on) . Adopt- 
ing their model of an inclined disc of enhanced wind, the 
Gorenstein (1975) approximation suggests the maximum X- 
ray column observed is consistent with a maximum value 
of M ~ 1CT 5 M Q yr _1 along the line of sight to the X-ray 
source. If the disc-wind has a solid angle of 10 per cent of 
An then, to be consistent with the ISO data and maximum 
column, a mass-loss rate of 1 x 1CP 3 M© yr _1 sterad -1 in 
the disc-wind is required; for 25 per cent of 47r it is 4 x 1CF 4 
M Q yr _1 sterad -1 ; and for 50 per cent of 4-zr it is 2 x 10 -4 
M Q yr -1 sterad -1 with 10 -5 M Q yr -1 outside the disc; i.e. 
disc density enhancements of 100, 40 and 20 respectively. 

We compare these results with previously reported the- 
oretical and observational work on deviations from spherical 
symmetry and density enhancements. The pioneering work 
of Schmid-Burgk (1982) demonstrated that deviations from 
spherical symmetry did not change the spectral slope but led 
to over-estimates of the mass loss rate. However, the maxi- 
mum correction factor found was only a factor of 2, for a pan- 
cake (disc) with a 10 to 1 aspect ratio. Ignace, Cassinelli & 
Bjorkman (1996) investigated density enhancements caused 
by rotation. For a WN5 wind their model found a density 
approximately 3 times greater than the global average in an 
equatorial cone of 10 degrees half-angle (or pole-to-equator 
density ratio of about 6). Rotation at only 10 per cent of 
critical was required. 

Harries, Hillier & Howarth (1998) found that 15-20 per 
cent of WR stars have significant polarization, which they 
attribute to global wind asymmetries, found only in the 
fastest rotators. They also found that these stars have the 
highest values of M. They postulate that these stars have 
equator to pole density ratios of 2 or 3 (with M scaling 
linearly with density). St Louis et al. (1995) propose a sec- 
tored wind structure for the WN star WR6. They find the 
wind divides roughly 50:50 into fast and slow wind regions, 
with the higher continuum fluxes coming from the hotter, 
higher speed regions. Radio imaging of the WN8 star WR147 
(Williams et al., 1997) report a wind 'shape' with a 3 to 2 
aspect ratio. There have been suggestions of a disc in the 
WR6 star WR134 - see Vreux et al. (1992), and references 
therein. Nota et al. (1995) report a bipolar outflow in the 
nebula around the WN8 star WR124, with an opening angle 
between 30 and 45 degrees. The binary nature of Cyg X-3 is 
likely to create deviations from spherical symmetry, as pro- 
posed originally by van Kerkwijk (1993) and for WR140 by 
Williams et al. (1990). The disc configuration reported by 
White & Becker (1995) is disputed by Harries et al. (1998). 



If the WR star's rotation is tidally locked to the 4.8 h 
orbital period, then it will be rotating at close to its criti- 
cal velocity and an equatorial enhancement of the mass loss 
rate is likely. We conclude that deviations of the wind from 
spherical symmetry need to be taken into consideration. Ig- 
noring such asymmetry tends to overestimate the the mass- 
loss rate. However, to be consistent with the visibility of 
the X-ray source the density contrast required by our re- 
sult appears to be large, and we are not aware of any work 
quantifying the overestimate for such large disc density en- 
hancements. 

4-1.2 M from P - an examination of assumptions 

To obtain the total stellar mass-loss rate which drives the 
orbital changes any other mass-loss rates occurring need to 
be added to this value of wind mass-loss rate. The X-ray 
luminosity of Cyg X-3 implies an accretion rate on to the 
compact object of ~ 10 -8 M© yr -1 , which is small compared 
with the mass lost through the wind. Mass will also be lost 
into the radio jets, but this is an uncertain quantity, and the 
mass in the jets could be anything between 0.1 to 100 per 
cent of the in-flowing mass (Fender & Pooley 1998, 2000). 
These additional masses are in the sense to exacerbate the 
problem. 

In obtaining M from P or b several other assumptions 
are made: that there is no loss of angular momentum other 
than the specific orbital angular momentum of the WR star 
carried away by the wind; that the rotating bodies are rigid; 
that the 4.8-h modulation, whatever it is, is rigidly locked to 
the orbital period; that the orbital period is 4.8-h, and not 
some multiple of that figure; that the eccentricity is zero; and 
that there are no other dynamical effects (e.g. gravitational 
radiation) or torques (e.g. magnetic) affecting the value of 

Porb • 

The assumption that the eccentricity is small seems jus- 
tified - given the amount of circumstellar matter a circular 
orbit seems likely, and no change in the light curve which 
could be attributed to precession has been observed. For 
convenience we adopt e = 0. It seems likely that 4.8 h is 
the fundamental period, but until its nature is understood 
there remains the possibility that the 4.8 h modulation is not 
rigidly tied to the orbital period and that the P observed is 
not accurately that of the binary system. 

Additional angular momentum could be lost through 
the jets, or through magnetic torques, for example, and ig- 
noring such losses would cause P or b to be underestimated. 

Gravitation radiation results in the orbital size shrink- 
ing, so neglect of this effect will cause the P or b due to the 
mass loss to be underestimated. The effect is most extreme 
for the fastest orbital velocities which occur for the most 
massive components (assuming the eccentricity is zero). 
Taking a Wolf-Rayet mass of 10 M and a black hole mass 
of 40 M , Pgr ~ -6x 10" 11 , or about 10 per cent of the 
observed P. If the orbital eccentricity were 0.5, the effect 
would be approximately 5 times larger, so although gravita- 
tional radiation produces an effect with the correct sign, it 
is not probably large enough to explain the discrepancy. Or- 
bital motion in a viscous medium (the wind) will also cause 
the orbit to shrink. 

It is known that tidal effects in early-type stars give rise 
to gravity waves which are radiatively damped, leading to 
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circularisation of elliptical orbits and the synchronisation of 
orbital and stellar rotation rates (Kumar & Qutaert 1998). 
Kumar and Qutaert show that if the stellar rotation rate is 
less than the orbital rotation rate (strictly periastron rate), 
P or b is negative. If we assume the Wolf-Rayet star rotates 
like its predecessor O star, then its rotational period would 
be a factor 10-20 times longer than the (current) orbital 
period. Tidal forces will presumably have spun it up; if it is 
now (close to) synchronised then for a 3 Rq star it will have 
a mean equatorial velocity of ~ 750 km s _1 and, assuming 
a mass of 10 Mq, be rotating at 95 per cent of its critical 
velocity. It seems likely that there have been tidal forces 
which have tended to contract the orbit and diminish the 
observed P or b, but in the absence of information about the 
current rotational period of the Wolf-Rayet we are unable 
to estimate the current magnitude of this effect. 

In summary, there are several very probable effects 
which could cause P or b to be underestimated - namely un- 
recognised sources of angular momentum loss, tidal effects, 
and to a lesser extent energy loss through viscous drag and 
gravitational radiation. Of the two methods of determining 
the mass loss rate, this seems to have the larger number of 
uncertainties. 



4.2 Flare data 

It is clear from the observations in 1997 June (Fig. | & 
^ that the ISO data reveal bright flaring at 4.5 and 11.5 
/xm at the same time as enhanced radio and X-ray activ- 
ity (while not presented here, the X-ray monitoring data, 
which reveal an X-ray bright state at this time, are pre- 
sented in FHP99). In addition, the spectrum appears to be- 
come flatter (although we caution that the lack of absolute 
simultaneity of the observations in the different filters adds 
additional uncertainties to the interpretation of spectral in- 
dices). This flattening is in qualitative agreement with the 
observations of Fender et al. (1996, 2000) who find that dur- 
ing both rapid (minutes timescales) and more gradual (days) 
variability, Cyg X-3 generally becomes redder in the infrared 
as it brightens. 

While Fender et al. (1996) interpret the flaring compo- 
nent from Cyg X-3 in terms of thermal free-free emission 
from a hot (T > 10 6 K) gas, infrared oscillations from GRS 
1915+105 (Fender & Pooley 1998 and references therein) 
have been interpreted as being synchrotron in origin. Either 
way, the rapid, large-amplitude, broad-band variability in 
the infrared fluxes during outburst periods implies the oc- 
currence of extremely violent processes on relatively small 
physical scales. While the timescales do not preclude a sim- 
ple change in the brightness of the thermal (wind) emis- 
sion, the significant spectral flattening appears to rule this 
out. Instead, similarly to GRS 1915+105, Cyg X-3 appears 
to display an approximately flat-spectrum (a ~ 0) com- 
ponent during outbursts. This may be related to the flat 
spectral component now also known to be present in Cyg 
X-l (Fender et al. 2000). The nature of this component is 
uncertain, but it appears to be related to jet formation, and 
possibly to instabilities in the accretion flow/disc, at least in 
GRS 1915+105. In summary, while there may be variability 
in the strength of the underlying thermal (wind) compo- 
nent, we believe the additional flux observed in outburst to 



be dominated by a variable, flat-spectrum component which 
is probably associated with the jets. 



5 CONCLUSIONS 

We report the first quiescent and flaring 4-16 fim observa- 
tions of Cygnus X-3. 

The quiescent spectrum is a continuation of the shorter 
wavelength IR spectrum, and fits well the WB75 model for 
thermal free-free emission from a wind. We obtain a mass 
loss rate, assuming spherical symmetry of ~ 10 -4 Mq yr _1 , 
and note that while this is consistent with other determi- 
nations, it is an order of magnitude larger than both the 
mass loss rate required to derive the change in orbital pe- 
riod (assuming Mr of the order 10 Mq), and the mass loss 
rate expected from the hydrogen column density to the X- 
ray source. These measurements are based on assumptions, 
which have been scrutinised carefully; while we can identify 
a number of areas where factors of a few can be gained, 
we have been unable to identify a single area where a fac- 
tor of ten can be found. We agree with the conclusion of 
FHP99 that the wind is most likely spherically asymmet- 
ric, but if this is the only effect operating, large disc density 
enhancement factors are required. There are several phe- 
nomena which could cause P or b to be underestimated; we 
regard this determination of M to be less reliable. 

During radio and X-ray outburst, the broadband in- 
frared spectrum shows rapid variability and a spectral flat- 
tening associated with increased flux densities. This is in 
agreement with what is seen in the near-IR during rapid 
(minutes) flaring events (Fender et al. 1996). Similar (but 
not identical) near-IR events observed from GRS 1915+105 
appear to be associated with jet formation episodes. The 
emission mechanism responsible for this flat spectral com- 
ponent remains uncertain; both free-free and synchrotron 
emission have been proposed already in the literature, and 
await detailed studies to test their viability. Whatever the 
origin of this component, it seems likely that it is associated 
with hot gas in the jet from the system, which naturally ex- 
plains its emergence during periods of enhanced radio (jet) 
activity. 
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